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ERGODICITY OF TWO DIMENSIONAL TURBULENCE 

ANTTI KUPIAINEN 



Abstract. We review the work of Hairer and Mattingly on ergodicity of two dimen- 
sional Navier-Stokes dynamics and discuss some open mathematical problems in the 
theory of 2d turbulence. 



1. Introduction 



The problem of turbulence has been described as the last great unsolved problem of 
j^ . classical physics. Understanding of the complicated motion of fluids in the presence of 

obstacles or stirring has been a challenge to mathematicians, physicists and engineers 
for quite a time now. The equations governing macroscopic fluid motion, the Navier 
Stokes equations, have been known for close to two centuries. For an incompressible 
fluid in units where the density equals one they read 



(1) dtu + u ■ Vu = z/Am — Vp + /. 

u{t,x) G Mf^ is the velocity field at time t at a; G A, a domain in Mf^ subject to the 
incompressibility condition 

(2) V ■ n = 

^ ■ and suitable boundary conditions on dA. v is the viscosity coefficient of the fluid, 

t^^ , p{t,x) the pressure and f{t,x) the external force that sustains the flow. Given / and 

^(0, ■) the task is to flnd u and p. It is fair to say that theoretical understanding of 
O ' the consequences of these equations is still in its infancy. On the mathematical side, 

l/^ ■ existence of smooth solutions for the three dimensional NS equations is wide open and 

O I has been chosen by some as one of the major problems of mathematics |Clay| . On the 

physical side, experimental violations of the Kolmogorov scaling theory of turbulence 
|Fri] are still waiting for theoretical understanding. 
•'-j . In two dimensions, i.e. for flows on the plane, there has been some progress during 

i^ ! the last ten years. On the physical side, 2d turbulence has been the subject of accurate 

c^ I numerical and experimental studies |Bof] . |RDCE] and mathematically the ergodic 

theory of the NS flow has been under intensive study. 

It is important to realize that for the problem of turbulence one is interested in a 
very particular kind of force in ([1]), namely one that has a flxed length scale L built 
into it. Examples of this are flows past obstacles, with L the characteristic size of 
the obstacle. In such a setup the flow exhibits universal statistical properties as the 
viscosity parameter tends to zero (actually the control parameter is a dimensionless 
quantity, the Reynolds number given by — where u is a velocity scale related to the 
forcing). E.g. time averages of measurements of suitable functions of u seem to show 
statistical properties only depending on the Reynolds number. It is therefore of some 

interest to inquire about the foundations for such statistical studies i.e. about the 
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ergodic properties of the NS flow in the turbulent setup of a fixed scale high Reynolds 
number forcing. 

A convenient model for isotropic and homogenous turbulence (i.e. in the limit of 
large Reynolds number and away from the boundary dA) is to consider □ eq. ([1]) on 
the torus T^ = M^/(27rZ)^ and take / random, a Fourier series with finite number of 
terms and coefficients independent white noises (see below). Then the deterministic 
dynamics of ([T]) is replaced by a Markov process and one may pose questions on its 
ergodic properties: whether the process has a unique stationary state and whether this 
is reached and with what rate from arbitrary initial conditions. 

This Markov process is a diffusion process of a very degenerate type. While the phase 
space is infinite dimensional the noise is finite dimensional. There are two general 
mechanisms that can contribute to the ergodic and mixing properties of stochastic 
fiows. One is dissipation, coming in our case from the Laplacean in ([1]). Dissipation 
contributes to ergodicity by exponential contraction of phase space under the fiow. 
Second mechanism comes from the spreading of the noise from its finite dimensional 
subspace due to the nonlinear term in ([1]). In finite dimensional diffusion processes 
this leads to hypoellipticity if the noise spreads to the full phase space: the transition 
kernels are smooth (for equations with smooth coefficients). Combined with some 
irreducibility of the process ergodicity follows. 

In our infinite dimensional setup the dissipation due to the Laplacean leads to strong 
damping of large enough (depending on the Reynolds number) Fourier modes. If we 
keep noise on all the other, low, modes then one can reduce the problem to a low 
mode dynamics, albeit with some (exponentially decaying) memory due to the large 
modes. Proofs of ergodicity and mixing of the dynamics were given in this case in 
the works |BKL] . |EMS] and |KS] . However, it seemed far from trivial to extend the 
hypoellipticity ideas to the infinite dimensional setup to control also the case of very 
degenerate forcing where the number of forced modes does not depend on the Reynolds 
number. This was accomplished by Hairer and Mattingly |HM1] . |HM2] who gave sharp 
sufficient conditions for the noise to produce ergodic and mixing dynamics. In what 
follows I will present the main points of their approach focusing on the difference to 
finite dimensional hypoelliptic diffusions. The papers |HM1] . |HM2] are very clearly 
written and they contain plenty of background material, especially [HM2j which also 
builds a more general formalism applicable also to some react ion- diffusion equations. 
|HM2] also corrects a mistake in |HM1] so it should be consulted for a thorough study. 
In the final section I discuss more informally what we have learned about 2d turbulence 
and what issues might be accessible to a rigorous mathematical analysis. 

I would like to thank J. Bricmont, M. Hairer and J. Mattingly for comments on this 
exposition and the European Research Council and Academy of Finland for financial 
support. 

2. 2D NS EQUATIONS 

The fundamental fact that is behind both the mathematical and physical under- 
standing of 3d NS equations is energy conservation: in the absence of forces smooth 
inviscid fiow preserves the L^ norm oi u{t,-). In two dimensions there is a second 



To get to the turbulent state one actually has to modify ((T]) a bit, see Section 8 
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conserved quantity, the enstrophy, which is related to the H^ norm and which leads to 
quite different physics and to much better regularity. 
Let us first define the vorticity 

W = V X M, 

which in (i = 2 is a (pseudo) scalar: u = diU2 — d2Ui. The NS equation becomes in 
terms of a; a transport equation: 

(3) to = uAu — u ■ Vw + g. 

where g = ^1/2 — d2fi- We will assume the average force vanishes i.e. / f{t, x)dx = 0. 
Then ([1]) preserves the condition Ju{t,x)dx = which we will assume. The incom- 
pressibility condition ([2]) allows to write u = Au where the linear operator A is given 
in terms of the Fourier transform by 

(4) Aw{k) = i{k2,-h)k-^u{k) 

for keZ^\0. 

The enstrophy £ is defined to be (half of) the L^-norm of u: 



£ = 2 oj{t,xYdx := -||a; 



For a smooth u the condition V ■ m = leads to the absence of contribution from the 
nonlinear term to the evolution of the enstrophy: 

(5) — = —V I iVioYdx + / ugdx. 

where the first term on the RHS can be interpreted as an enstrophy dissipation rate and 
the second one as an enstrophy injection rate. Using Poincare inequality ||Va;|| > ||a;|| 
and simple estimates one deduces 

(6) ||u;(t)f <e-^*||u;(0)f + ^-^up||^(t)f. 

t 

This a'priori estimate for the H^ norm of u is the main ingredient in the proof of global 
regularity of the 2d NS flow. 

We wish now to discuss a version of ([3]) where the force g is random. We work in the 
subspace of real valued L^(T^) functions with a;(0) = 0. It will be convenient to use 
the following basis for this space. Let Z~^ be the "upper half plane" in Z^ consisting 
oi k = {ki,k2) with fcg > or fcs = and ki > 0. Hence Z^ \ = Z+ U {-Z+). 
Let Cfc = sin /ex for k G Z^ and e^ = coskx for —k G Z~^. For each A; G Z^ pick 
independent Brownian motions f3k(t) with unit speed, denoted collectively by /3(t) and 
numbers 7^ G M. Let 

(7) g/3(t) = $^7fc/3fe(t)efc. 

fcgZ2 

The stochastic version of eq. ([3]) reads 

(8) du = (z/Aw - u ■ Vu)dt + Qdl3. 
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Regularity of the stochastic flow proceeds in parallel with the deterministic case as long 
as 7fc have enough decay at infinity. The analog of the enstrophy conservation eq. ([5]) 
is obtained by an application of the Ito formula 

(9) d£ = ld\\ujf = -u\\Vujfdt + (tu, Qd(3) + edt 

where e = 27r^ ^^ 7^ can be interpreted as the enstrophy injection rate. Taking aver- 
ages we get a probabilistic analog of (|5]) and ([6]): 

(10) —E£ = -i^E\\Vuf + e. 

and 

(11) E||cj(t)f <e-2-*||u;(0)||2 + z/-ie. 

Actually ([9]) can be used to control exponential moments of the enstrophy |BKLj . |HM1] 
Lemma A.I.: 

(12) Eexp{r]\\uj{t)f) < 2exp(r7e-'^*||w(0)f ) 

for all ?7 < - i.e. probability for large L^ norm is exponentially small. ([9]) also allows 
to control the time integral of the H^ norm: 

(13) Eexp(r/z/ / \\Vu{t)f) <2exp{r]et + 7]\\u{0)f) 

Jo 

again for all r^ < -. Such a'priori estimates allow one to prove the existence and 
pathwise uniqueness of strong solutions to eq. ([8]) under quite general conditions on 
the noise coefficients 7^, see e.g. |Fla] and |MR] . Of course the PDE regularity is the 
harder the less decay the 7^ have at infinity. As explained in the introduction, for the 
turbulence problem only finite number of 7^ are nonzero. Thus from the point of view 
of regularity the turbulent case is easy (this is not true in 3d!). However, the ergodicity 
of the fiow becomes the harder the less there is noise. 

3. INVARIANT MEASURE 

Let us now specialize to the case where 7^ = for A; ^ ii' where K is a. finite 
set. Thus the noise is finite dimensional: /3 = {Pk}k£K can be identified with the 
Wiener process in fi = C([0, oo),M-^) where D = \K\, equipped with the Wiener 
measure W{db). The solution of eq. ([8]) is a one parameter family of continuous maps 
$i : fi X L2(T2) -^ L2(T2) such that uj{t) = ^t{f3,uJo) solves eq. (^ with initial 
condition uq and noise realization /3. Actually, $i is (Frechet) differentiable in /3 and 

Wo- 

u){t) is a Markov process with state space H = L^(T^). It gives rise to transition 
probabilities Pt{ujQ,A) which are probability measures on H, giving the probability of 
entering the set A G H at time t given that at time we have Ci;(0) = Uq: 

Ptiuo,A) = ElAiuit)). 

The transition probabilities generate a semigroup Vt on bounded measurable functions 
on H by the same formula: 

(14) Pt<P= [ Pti;du)<l){u) 
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and the adjoint semigroup acting on bounded (Borel) measures: 
(15) P:fi= f fiiduo)PtiuJo,-) 



We are interested in the invariant (or stationary) probabihty measures fi* satisfying 
the equation 

(16) P>* = /i*. 

Existence of an invariant measure is straightforward given the strong probabihstic 
control of the flow. One considers the family of time averages /i^ = t~^ /q dsPs{uJo, ■) 
and shows it is tight. Prohorov's theorem then yields a limit point which is shown to 
be invariant. 

Uniqueness of the invariant measure is much more subtle. It implies ergodicity, i.e. in 
particular the equivalence of time averages and ensemble averages: limt_j.oo /Wt " (0) = 
n{4>) for all G L'^{H, fi) and /i-a.s. in uq. In practice one would like to have more i.e. 
the convergence in some sense of the measures Pt{u;o, ■) to //* as t — )■ cxd. This leads to 
various mixing concepts. 

4. DISSIPATION AND SMOOTHING 

For finite dimensional diffusion processes it is well known that uniqueness of the 
invariant measure follows from recurrence and smoothing properties of the transition 
probabilities. Let us sketch a special version of this argument having the application 
to NS in mind. 

The semigroup Vt is called strong Feller if the image is continuous for (p measurable. 
This has drastic consequences for the supports of invariant measures. Recall that x 
belongs to the support of a finite Borel measure /i on a Polish space (our setup) if 
/i(t/) > for all open U containing x. Then the supports of two distinct ergodic 
invariant probability measures for a strong Feller semigroup are disjoint. To see this, 
suppose ^ -L u and x G supp/x fl suppz/. Pick A with fi{A) = 1 and i^{A) = 0. By 
strong Feller there exists a U containing x such that sup^^^g^ \Pt{y, ^) — Pt{z, A)\ < ^. 
Moreover by assumption a := min{/i(t/), i^{U)} > 0. Write /x = (1 — a) ft + a^jj with 
/2 and [lu probability measures with fJ.u{U) = 1 {i.e.fiu = fJ'^u/f'{U)) and u similarly. 
Then, by invariance \n{A) - i^{A)\ = IPX A) - Vtu{A)\ and thus 

1 = \f,{A)-u{A)\<{l-a)\V:mA)-V:HA)\+a\V:MA)-V:MA)\ 

< {l-a) + af \Vt{y,A)-Vt{z,A)\^iu{dy)Mdz)<l-\a 

JuxU 

a contradiction. 

Suppose now that we knew that there exists an x that necessarily belongs to the 
support of every invariant measure of a strong Feller semigroup. We could then con- 
clude uniqueness. This is a reasonable strategy for the NS equation. Indeed, u = 
is such a point. This follows since NS equation is dissipative. Without forcing the 
fluid slows down i.e. the L^ norm decays exponentially (see eq. ([6])). There is a non 
zero probability for the force to stay small enough so that any neighborhood of can 
be reached. More precisely, the Ito formula flTO|l combined with Poincare inequality 
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||u;|| < ||Vu;|| yields 

fi{du)\\u\\^ <e/u 

for every invariant probability measure fi. Hence, there exists R < oo such that every 
such measure has at least half its mass in the ball B^ of radius R centered at in H. 
Thus one needs to show: for all r > there exits T^ < oo such that 

Ir:= inf Pt^{coo, Br) > 

(see HM], Lemma 3.1.). Then /i(5^) = PX^^) > Ur > for all r > 0. 

This strategy does not quite work in our case since the strong Feller property is very 
hard to show for Vt and might very well not be true. One of the main accomplishments 
of Hairer and Mattingly was to replace it with a condition that is more natural for NS 
and yet allows one to conclude that the supports of invariant measures are disjoint. 

5. ASYMPTOTIC STRONG FELLER PROPERTY 

A strong Feller semigroup maps bounded functions to continuous ones. Often the 
easiest way to prove this is to show a bit more |DPZ] . Lemma 7.1.5.: 

Proposition 5.1. A semigroup on a Hilbert space H is strong Feller if for all (p : if — t- M 

with \\4>\\oc '■= ^'^PxeH \^i^)\ ^^<^ ||-D0||oo finite one has 

(17) \\DVt<P{x)\\<C{\\x\\)UU, 
where C : M^ — t- R and D is the Frechet derivative. 

We will now argue that the condition (TT7I) is not very natural for the NS dynamics. 
As mentioned in the introduction there are (at least) two ways ergodicity can result. 
One is due to smoothing by the noise, the other is due to dissipation that erases memory 
of the initial conditions. The former effect leads to a condition like ( TT7|1 . the latter not. 
Let us next discuss the latter effect in our case. 

Let Js^t with s < t he the derivative of the NS flow ([8]) between times s and t, i.e. 
for every C, E H, Jg^tC, '■= (,(t) is the solution of the linear equation 

(18) dt^t) = i^^^it) + Aujit) ■ Ve(t) + A^t) ■ Va;(t) := C^it)m 

ioT t > s and ^(s) = ^. This linear equation is readily controlled in terms of the H^ 
norm of a; ( [HMlj . Lemma 4.10): 

(19) \\^(t)\\<exp{C{6,u){t-s) + 6 f \\Vuj{r)\\'dr)U{s)\\ 

J s 

for any 6 > 0. Combining with the a priori estimate (TT3l) then 

(20) nJs,tV < 2^exp (C(e, iy,r],p){t - s) + vMs)f) 

for all ?7 > 0, all p < oo. 

Eqs. (fT9l) and (!20|) indicate possible exponential separation of trajectories. However, 
since the Laplacean is the Fourier multiplier —k"^ it is not surprising that the high 
Fourier modes of ^ are strongly damped for a time that can be taken as large as we 
wish as A^ is increased. This is expressed by [HMlj . Lemma 4.17.: 
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Lemma 5.2. For every p > 1, every T > 0, and every two constants 7, ?7 > 0, there 
exists an orthogonal projector vr^ onto a finite number of Fourier modes such that 

E(||(l - vr,)Vr + II V(l -7r,)r) < 76^'"'"°"' 

For such contracting dynamics fITTj) is not a natural condition to try to prove. Indeed, 
let ^h = {^~TTe)C, be the projection of ^ to the high modes and consider the toy problem 
where we apply (1 — iri) to eq. flTSj) and drop altogether the w-dependent terms: 

dMt) = i^^Ut)- 

Then for a function 0(ci;) = '?/'((l — 7ii)u}) depending only on the high modes we have 
DVt(p{coo)^ = ED^{u{t))^hit). Since in this toy case ||^h(t)|| < e-^*||^|| for A > we 
conclude 

\\DVt(p{x)\\ < e-^'\\D<f)\U 

This toy model and Lemma 5.2. motivate the following definition by Hairer and Mat- 
tingly ( jHMl] . Proposition 3.12.) 

Definition 5.3. semigroup Vt on a Hilbert space "H is asymptotically strong Feller 

if, there exist two positive sequences in and 6n with {t„} nondecreasing and {5„} 
converging to zero such that for aii : "H — t- M with \\4>\\oo SLnd \\D(f)\\oo finite, 

(21) \Dn„cP{x)\<C{\\x\\){U\\^ + 5„]m\\^) 
for all n, where C : R+ — ;• M. 

(Hairer and Mattingly actually give a "topological" definition of asymptotically strong 
Feller condition which is implied by the one above). The main point is the following 
result whose proof is similar to the one given above in the strong Feller case ( [HMlj , 
Theorem 3.16): 

Proposition 5.4. f the semigroup is asymptotically strong Feller at x then x belongs 
to the support of at most one ergodic invariant measure . 

We saw above that it is not unreasonable to expect that the high mode dynamics 
give rise to the second term in fl2T]) . Thus the question remains: why would the low 
mode dynamics be strong Feller? The answer to this question lies in the hypoellipticity 
of the low mode dynamics. 

6. HYPOELLIPTICITY 

Let us think about the low mode dynamics first in the Galerkin approximation i.e. 
by putting the high modes to zero. More formally, consider the equation 

(22) du = {uAu - TTeiu ■ Vuj))dt + Qdp. 

where we assume the forcing is on low modes (1 — 'n:e)QP = and set (1 — 7i£)u) = 0. 
Eq. (I22!) defines a diffusion process in a finite dimensional space which we may identify 
with M^, A^ = dimiTgH. The diffusion process is thus degenerate with the dimension 
D of the noise (much) smaller than N. The strong Feller property follows for such 
diffusions provided the generator of the diffusion process is hypoelliptic. Let us discuss 
this next. 
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Recall the Fourier basis {ck} for H. Let the range of vr^ be the span of {e^} with 
\k\ < M. Write u = ^^Wfcefc. Then the equation fl22|) reads 

(23) duk = Vk{uj)dt + -fkd(3k. 

where Vk is given by 

j+e=k ' ' '"'' 

and Wk = h^-k + ^i^k for k G Z^ and w^k = ^k- The generator of this diffusion is 



(24) L = Xo + J2^. 



k&K 



where we recall that 7fc = for k ^ K. The vector fields Xa are given by 

^0 = ^ Vkdui^ 

k 

Xk = ikd^^ 

An operator of the form fl24|) with smooth vector fields X^ is known to generate a 
semigroup Vt with smooth kernel (hence it is strong Feller) provided the Hormander 
bracket condition is satisfied (L is then hypoelliptic). The condition is that the span of 
the vector fields X^ , j 7^ and [Xi^ , [Xi^ , . . . [Xj^.^ , Xj^.]] . . .] ioi k > 1 and ij G {0} U K 
at each cu G M.^ equals R^ . 

To check this condition in the NS case is a purely algebraic excercise and the result 
is the following jEMllHMl] : 

Proposition 6.1. he following conditions for the set K ^1?\ {0} are sufficient for 
the Hormander bracket condition to be satisfied: 

(a) K is invariant under the reflection k -^ —k 

(b) K contains at least two elements of unequal length 

(c) K spans T? under linear combinations with integer coefficients. 

An example of a very degenerate forcing that suffices is given by the the set K = 
{(1, 0), (—1, 0), (1, 1), (—1, —1)} i.e there is forcing only on two wave vectors and their 
reflections. 

Note that the Proposition 5.1. is true for arbitrary (large enough) Galerkin cutoff N. 
Hence the full infinite dimensional generator formally satisfies the Hormander condition 
and one might be tempted to try to use this to return to the attempt to prove the strong 
Feller property for Vt- However, it is likely that, as we let N increase, the derivatives 
of the kernel of Vt with respect to the high modes blow up since the smoothing is very 
weak for them. It is much more natural to try to use in that regime the dissipation as 
coded in the asymptotic strong Feller condition. 

Let us finally remark that if all the 7^ in (l23l) are nonzero the generator L is elliptic. 
If N is large enough (of the order e/u^) then one may use the dissipativity of the high 
mode dynamics to solve for the high modes in terms of the (temporal history) of the 
low modes and use the ellipticity of the latter to prove ergodicity and mixing of the 
full dynamics pm lEMSl iKS] . 
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7. MALLIAVIN MATRIX 

Why does elliptic diffusion produce smoothness in transition k;ernels? One way to 
think; about this is to consider trajectories of the flow. Noise will make the trajectories 
non-unique: a change in the initial condition can be compensated by the noise. In 
elliptic diffusions noise spans the whole space and the compensation is immediate, in 
hypoelliptic diffusions the nonlinearity spreads the noise in all directions thanks to the 
bracket condition. Thus a derivative of the solution in the initial condition should 
equal its derivative in a particular direction in the (history of) noise space. Since we 
are integrating over the noise the latter derivative can be integrated by parts and hence 
an estimate like the strong Feller property can emerge. 

To be more explicit recall that we wrote the solution of the stochastic NS equa- 
tion as u(t) = ^t{f3,o-'o) with <I>t smooth in the noise (3 G C([0, oo)),]R-^) and the 
initial condition cjq- Also we have denoted the derivative in the initial condition by 
{D^Mt),0 = Jo4 = ^W for ^eH. Thus 

(25) {DVt<Piuo),0 = E((Z}0)(a;(t)),e(t)). 

Consider next the infinitesimal change in the solution corresponding to the change of 
the noise /3 in the direction V e C([0, oo)),R^): {Dpu{t), V) := ({t). ({t) satisfies the 
same linearized NS equation flTSj) but with with forcing QV: 

(26) dC{t) = C^it)C{t)dt + QdV{t) 

and zero initial condition. The natural space to vary the noise is the Cameron-Martin 
space i.e. to take V of the form V{t) = /^ t'(s) c/s with v G ^^^^([0, oo],M^). By 
variation of constants ({t) is then given by 



(27) C{t) = / Js,tQv{s)ds := Atv. 

Jo 

Actually, the v one will eventually use is itself a function of the noise (see eq. (l40|) ). 
but it will be a.s. in Lf^^. The upshot is that At : -^^([0, t], M^) — )■ if is an a.s. bounded 
random operator and so the Frechet derivative can be written as 






fee/c 



where the operator V^ is called the Malliavin derivative and heuristically corresponds 
to an instantaneous kick at time s to the direction k in noise space. Explicitly 

(28) V':co{t) = Js,tlkek. 
Suppose now we can find a v such that 

(29) m = C{t) i.e. Jo,t^ = Atv 
Inserting this to eq. fl25l) we get 

(30) {DVt(l){uo),0 =mD4>)i^it)),at)) = nD^^{co{t)),V). 

The derivative Dj^ in eq. (!30l) can be integrated by parts in the Gaussian Wiener 
measure to obtain 

(31) E{Dp^{u{t)), V) = E{(P{co{t))D;V). 
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In other words, the expression D^ is the adjoint of -D^ in L^{Q, W). If the process v is 
adapted to the Brownian filtration its expression is simply D*^V = J2k I ''^k{s)d/3k{s), 
the Ito integral. Otherwise a derivative of v with respect to the noise also appears and 
DIV is called the Skorokhod integral off. Combining ( 13T]) with (!30|) the desired bound 
follows: 

(32) \{DVt<P{uo),0\<\moonD;V\. 

It remains to solve eq. fl29l) for V (i.e. for v). Let A^ be the Hilbert space adjoint of 
At i.e. explicitly 

(33) {a;o{s) = q*ju 

for s < t. Then the Malliavin matrix is defined by 

(34) M{t) := AA; = [ Js,tQQ*Jltds 

Jo 

Suppose M(t) is invertible. Then, clearly a solution to (129|) is given by 

(35) v = a;m^JJo4 

To sketch the rest of the story in the finite dimensional setup we need a bound for 
the Skorokhod integral appearing in eq. (132|) [Nua] : 

(36) E{D*pVf <E f |i;(s)pcis + Ve f V'lvi{r)V[vk{s))dsdr. 

"^0 kl ^ 

The first term is the usual identity for the L^ norm of the Ito integral, the second term 
appears for a non-adapted v, as is the one given by (135|) . To compute the Malliavin 
derivative of v in (136|) note that all we need is to compute V^Js^t since At and M{t) 
are expressed in terms of Jg^t- This in turn is obtained by differentiating the equation 
(|T8|): r] ■.= V'^^ satisfies 

dtrj{t) = uA^it) + Auj{t) ■ Ve(t) + A^t ■ Vu;(t) := C^it)r]{t) + B{Jr,tlkek,v{t)) 

where B is the bilinear form in appearing in C. By variation of constants an expression 
involving only J emerges. 

Thus in the finite dimensional setup (so e.g. for the Galerkin NS) the main work 
to be done is to show that M{t)~^ has good probabilistic bounds. Indeed it turns out 
||M(t)~^|| is in L'P{Q) for all p < oo. In the infinite dimensional case with degenerate 
noise it is unlikely that M{t) is a.s. invertible. QQ* is proportional to the projection 
in H to the subspace generated by the noise. In the expression for M{t) the dynamics 
spreads the range beyond this subspace, however we expect the projection of the result 
to the high modes to be very small. The key estimate on the Malliavin matrix Hairer 
and Mattingly prove is that M{t) is unlikely to be small on vectors that have large 
projection to low modes: 

Proposition 7.1. or every a, t], p and every orthogonal projection 7T£ on a finite number 
of Fourier modes, tliere exists C such that 

(37) P( inf ^^^<e)<Ce^exp(r^||a;of) 



holds for every e G (0, 1), and for every uq G H. 
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We will not discuss the details of the proof which is the technical core of the paper 
|HM2] (see also |MP] ). However, for experts we want make the following comments. 
One major difficulty H&M face is that the integrand in the expression for the Malliavin 
matrix is not adapted, i.e. depends on the future noise. The usual way out of this 
problem in the finite dimensional theory is to use the semigroup property Jo,t = Js,tJo,s 
to rewrite M{t) = Jo^tM{t)J^t with 



Mit) := / Jo,lQQ*Jo,s' 
Jo 



the reduced Malliavin matrix (and the control v{s) = Q*jQg^^M(t)^^$,). In finite 
dimensions Jo,t is invertible and now the integrand is adapted. The proof then uses 
Norris lemma |Norj which states that if a semimartingale is small then both its bounded 
variation part and local martingale part are small. In the infinite dimensional case with 
degenerate noise Jo,i is not invertible due to dissipation of the high modes. Hence one 
needs to work with non-adapted processes. The way out for H&M is the polynomial 
nature of the nonlinearity. In the iterative proof to show that (</>, M{t)(j)) small implies 
s — )■ {Js^tP{u{s))(f), (f)) is small for the various multiple commutators P the P will always 
be a polynomial. One then writes u{s) = v{s) + QP{s) where v is more regular and 
expands P{u{s)) in powers of QP{s), ending up with a polynomial in the Wiener process 
/3(s) with coefficients that are nonadapted processes, but with higher regularity. The 
basic lemma one now needs is that such a Wiener polynomial can be small only if all 
the coefficients are small (up to events of small probability). 

8. LOW MODE CONTROL 

The approach to prove smoothness sketched in the previous section is a form of 
stochastic control where the noise is used to force solution to a prescribed region in 
phase space (for results on stochastic control in our setup see also [ASj . [AKSSj ). 
We saw that an exact compensation of the change of initial condition by a change 
in the noise seems impossible, but Proposition 6.1. gives reason to hope that partial 
compensation is possible for the low modes. Since by Lemma 5.2. the high modes are 
contracted the idea of Hairer and Mattingly is to do an approximate control such that 
instead of the full control fl29l) we have C,(t) — ((t) — )■ as t — t- oo. Thus, as before 
let V G L^Q^(R+,M^) be a shift in the noise and ((t) = v4o,if be the corresponding 
(infinitesimal) shift in the solution. Let 

(38) Pit) = at) -at). 

Then, instead of the identities fl30|) and fl3Tl) we obtain 
(39) {DVMuo),0 = E{^{u{t))D;V)+E{D(P{co{t)),p{t)) 

< U\UE\D;V\ + \\D<j)\\^E\\p{t)\\. 

The asymptotic strong Feller property will follow provided v can be chosen such that 
E|D^V^| stays bounded as t — )■ oo and E||p(t)|| tends to zero exponentially. 

To find V H&M use use a construction where at successive time intervals two steps are 
alternated, one where high modes contract, the second where low modes are controlled 
by the noise. Suppose at some time t we knew p(t) is mostly in the high mode subspace, 
i.e. ||7r£p(t)|| << ||p(t)||. Then, at least for a short time it pays to set f = since the 
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linearized dynamics contracts such a p strongly. However, we cannot do this for too 
long since the low mode part of p will increase. Then provided we can find a v that will 
compensate the low mode part on a fixed time interval while leaving the high mode 
part approximately intact we can iterate the procedure. 

The low mode control is a simple modification of the full control explained in the 
previous section. Let us take the time intervals as [ra, n + 1] with n odd integer for the 
first step and even integer for the second step. Thus we set v{t) = for t G [n,?! + 1], 
n odd. Let An := v4„_„+i, M„ := v4„y4* and J„ := J„,,„+i. For n even take 

(40) Vn := V\ln,n+l] = Al{Mn + f3)-^JnP{n). 

Note that except for the parameter /3 this agrees with the full control (!35|) . While for 
/3 = the inverse in fHOl) most likely doesn't exist, for /3 > it does. The point now is 
that for small enough /3 fHOl) does a good job for the low mode control while the high 
modes remain approximately intact. To see this, compute 

(41) p{n + l) = ^{n + 1) - C{n + 1) 

= Jni{n) - (JnCin) + AnA:^{Mn + /3)-V„p(n)) 
= /3{Mn + /3y'Jnp{n) 

By Proposition 6.1. eigenvectors of M„ with small eigenvalues have small projections 
to the low modes. Hence one expects that for small /3 the operator /3(M„ + /3)~^ is small 
on vectors ip with ||7r£'j/'|| > a||'?/'|| whereas it is obviously bounded by one elsewhere. 
Combining the two steps we get the iteration 

(42) p{n + 2) = Jn+i/3iMn + P)'^JnP{n) 

Combining Lemma 5.2, the bound fl20|) and Proposition 6.1. H&M prove ( [HMlj . 
Lemma 4.16) 

Proposition 8.1. or every two constants '^,1] > and every p > 1, there exists a 
constant /3o > such that for n even 

holds almost surely whenever /3 < f3o. 

Iterating Proposition 8.1. the exponential decay of E||p(t)|| then follows ( |HMlj . 
Lemma 4.13). 

What remains is to bound the term E|Do(5^| in fl39|) uniformly in t, i.e. to bound 
the two integrals in fl36|) . The crux of the matter here is that both terms can be written 
as a sum over n of factors proportional to p{n) which provides a convergence factor. 
For the first term this is obvious by (HOj) . For the second one we need to go back to 
the integration by parts formula eq. fl3T|) . By construction V(t) = f^ v{s)ds = J2n^n 
where Vn is J-'n+2 measurable. Thus since the integration by parts is local in time 
D*aV = V^ D% Vn and the second factor becomes 

(43) 5Z^ / ii{Vsv{r),Vrv{s))dsdr. 

For details of how to finish the argument we refer the reader to Section 4.8. in [HMlj . 
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9. TURBULENCE 

We have seen that the NS dynamics has a unique stationary state under very general 
conditions of the forcing. Moreover, it can be proven that the dynamics is mixing |HM3j 
and the stationary state is reached exponentially fast from arbitrary initial conditions 
and for arbitrary large Reynolds numbers R (for earlier proofs of mixing in the case 
where i?-dependent number of modes are forced see |BKLj . [Mat] ) . Does this mean we 
have reached some understanding on the properties of this state, in particular on the 
phenomenon of turbulence? The proof outlined in the previous sections uses properties 
of the system that have counterparts in the phenomenological theory of turbulence. 
These are the dissipation of the high Fourier modes and the transfer of the noise 
from the forced modes to the unforced ones due to nonlinearity. The latter point 
is significanta because most results of the NS dynamics are based on the energy and 
enstrophy conservation laws alone, and those bounds would hold even if the nonlinearity 
was zero. Therefore, the properties of the latter are not used. 

This being said it must be stressed that we have gained very little understanding 
of the actual nature of the invariant state. Crucial part of the proof is irreducibility 
which is based on the fact that a; = belongs to the support of every invariant measure. 
Recall that this holds, because there is a small probability that the random forces are 
close to zero for any given time interval so the fluid flow slows down due to viscosity. 
This is clearly not the true reason one sees fast approach to stationarity in physical 
experiments. The mixing times resulting from visits to the origin will be much larger 
than the ones observed. To understand the real mechanism for mixing one has to 
understand much better the transfer of energy and enstrophy from the forcing scale to 
other scales. 

It was Kraichnan's observation [Kra] that we should expect this transfer to be in 
two dimensions quite different from the three dimensional case. In three dimensions, 
according to the Richardson-Kolmogorov picture the forcing in low modes injects into 
the system energy which is transported due to the nonlinearity in NS equation to the 
higher modes and eventually dissipated by the viscous term by large enough modes. 
This transport of energy through scales in wave number space (i.e. \k\ := k) is called 
the Richardson energy cascade. In fact the theory predicts a constant flux of energy 
from the injection scale (in our case 1) to the dissipation scale Ky (these claims can be 
formulated in terms of various correlation functions in the putative stationary state, 
see e.g. the review [Kupj ). Kraichnan noted that the existence in 2d of the second 
conserved quantity of the inviscid flow, the enstrophy, means that one has to pose 
the question at what scales (if any) energy and enstrophy are dissipated and if there 
exist separate fluxes for the two. His observation was that the fluxes of energy and 
enstrophy are to opposite directions, energy flows towards low modes and enstrophy 
towards high ones. Moreover, energy tends to be not dissipated at all whereas enstrophy 
is dissipated at high modes like energy in the 3d case. The presence of the two cascades, 
the direct cascade of enstrophy and the inverse cascade of energy is very well established 
both numerically [Bof] and experimentally |RDCE] . In what follows we will point out 
a couple of mathematical questions regarding this picture which would be nice to 
understand. 

To state the Kraichnan picture more precisely it is convenient to work on a torus of 
size N i.e. T^ := (R/(27riVZ)^ rather than A^ = 1 we had before. Of course by simple 
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scaling we can get rid of the A^ at the expense of changing v and the forcing scale, 
but since the theory involves large separations of the scales of dissipation, forcing and 
injection it is natural to take N large (eventually to infinity) we rather not do that. 
Consider now the NS dynamics on T^ with the random forcing on Fourier modes of 
size \k\ ~ Kf » N~'^ (observe that now k G {N~'^ZY i.e. \k\ > l/N). We shall 
add to the NS equation ([3]) an extra term that damps the low Fourier modes more 
strongly than the viscous term does (note that uk"^ can be as small as u/N^). This 
is the Ekman friction term —tu for r > 0. Stationary states for this system exist for 
the same reasons as before and uniqueness should follow in the presence of the friction 
as without provided the conditions of Proposition 5.1. hold. The Kraichnan theory 
makes predictions on this stationary state, call it ij,u,t,n in the various limits A^ — t- oo, 
r -)■ and z/ -)■ 0. 

The starting point is conservation laws of energy and enstrophy following from the 
enstrophy balance equation (ITO|) and a corresponding one for energy and taking into 
account the extra friction term in the equation. Since the unique stationary state is 
translation invariant these become local identities, for enstrophy 

(44) z/E,,,,^(Vw(x))' + rE,,,,,v(w(x))2 = e 
and analogously for energy 

(45) z/E,,,,;v(Vm(x))2 + TE,^r,NH^)y = e' 

with e' the energy injection rate (per unit volume) which is proportional to e^^j^- Eu,t,n 
denotes expectation in the measure fiu,T,N- 

The first question to pose is what happens to the viscous dissipation of energy and 
enstrophy as z/ — )■ 0. All the evidence points to vanishing of energy dissipation 

(46) limz/E^^jv(VM(x))^ = 

Enstrophy dissipation is more subtle as we will see below, but again it is believed [Ber] 
that it vanishes: 

(47) lim uE^rNiVuix)? = 

It would be interesting to prove these statements and also to understand whether 
a limiting measure lim,^_j,o fJ'u,T,N exists and is supported on solutions of the damped 
randomly forced Euler equation. Indeed, some indications that this could be done 
comes from the work |CRj where time averages of solutions and statistical solutions 
are controlled in that limit. They are shown to be given in terms of solutions of the 
Euler equation and in particular |CRj prove the relation (1471) in that setup. 

The main predictions of the Kraichnan theory come from the limit N ^ oo and 
r — !• 0. The limit A^ — >■ oo means we are considering the NS dynamics in M^. It is an 
interesting problem to try to prove that the (weak) limit lim7v-s.oo fJ'u,T,N = f^u,T,oo exists. 
Note that we don't expect this state to be supported on L^ but rather on polynomially 
bounded (and presumably smooth) functions. The reason the large volume limit might 
exist is the damping of the low modes by the friction term. It produces an effective 
low wave number cutoff (which turns out to be ~ r^/^e'~^). 

Granting this, what happens if we now take r — )■ 0? Is there also a measure fiu,o,oo'^ 
The prediction of the Kraichnan theory is that the viscous energy dissipation P6|) van- 
ishes as z/ — 7- uniformly in r. Thus in that limit Ej,^t-.oo(m(x))^ = (e' — o{v))/t 
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i.e. the average energy density is not bounded in the putative hniiting measure 
^u,o,oo- However, it is beheved that /i^,o,oo is supported on smooth co and in partic- 
ular hmT-_^oTlEi,.T-,oo(i^(2;))^ = 0. Then (jUj) imphes dissipative anomaly ior enstrophy: 
enstrophy dissipation remains nonzero as i/ — )■ i.e. hm^^o ^^u,o,oo(^^{x))'^ = e > 0. 

The Kraichnan theory makes more quantitative predictions of the distribution of 
energy and enstrophy according to wave number. Define the energy spectrum for k G 



(48) e(/t) = 27ig{K)/K 

where ^(|fc|) is the Fourier transform of the vorticity 2-point function 

g{x -y) = ¥.y^r,oo(^{x)uj{y). 
Then energy density is given by 

nOO 



(49) E^,^,ooM(a;)^ = / e{K)dK 

Jo 

and enstrophy density by J e{K)n^dK.. Kraichnan theory predicts 



(50) e{K) ~ ^ 



g/2/3^-5/3^ K^ « K « 



K 



/• 



where n^, ~ z/^^es is the dissipation scale and k^ ~ x^/^e'~^ the friction scale. 

The picture painted by the Kraichnan theory on 2d turbulence is thus quite complex. 
With well separated scales of viscous dissipation, injection and friction energy flows 
from the injection scale towards small wave numbers and is eventually dissipated by 
the friction. In the absence of friction and in infinite volume energy flows to ever smaller 
wave numbers and energy density is not defined in the stationary state. Enstrophy in 
turn flows to high wave numbers and is dissipated there by the viscosity. Only in the 
state fj.u,o,oo as 1/ — )■ one expects to have constant fluxes of energy and enstrophy, for 
some exact calculations (subject to regularity assumptions) see [Ber] . One has to be 
careful with the order of limits as is seen from the behavior of enstrophy dissipation. 
Note in particular that the stationary state fiu,o,N which we have been discussing in 
the previous sections does not exhibit turbulence in the sense of cascades of energy 
and enstrophy. Here energy will reside in low modes, indeed, in experiments one 
often sees the formation of a few large vortices in the flow. If u is taken to zero in 
this state then both energy and enstrophy will blow up and indeed, no limit measure 
exists |Kukj . In [Kukj it is proven that only by taking the injection rate e (and thus 
also e') proportional to z^ a nontrivial limiting measure exists. Formally this limit still 
corresponds to diverging Reynolds number, but one does not expect it to be a turbulent 
state with near constant fluxes of energy and enstrophy. 

What makes the Kraichnan theory intriguing is that e.g. the spectrum (150|) seems to 
be very well verified numerically and experimentally. Moreover, the invariant measure 
seems to possess strong scale invariance properties, at least in the inverse cascade 
regime. There are even indications of conformal invariance |BBCFj . Thus it is not 
excluded that some of its properties could be mathematically accessible. 
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